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Abstract. In response to inconsistencies which have arisen in 
results from a hydrodynamic model in simulations of high ion 
temperatures (1-2 eV) observed in low density, outer plasmas- 
phere flux tubes, we postulate a reduced thermal conductivity co- 
efficient in which only particles in the loss cone of the quasi- 
collisionless plasma contribute to the thermal conduction. Other 
particles are assumed to magnetically mirror before they reach 
the topside ionosphere and therefore not to remove thermal en- 
ergy from the plasmasphere. This concept is used to formulate a 
mathematically simple, but physically limiting model for a 
modified thermal conductivity coefficient. When this modified 
coefficient is employed in the hydrodynamic model in a case 
study, the inconsistencies between simulation results and obser- 
vations are largely resolved. The high simulated ion tempera- 
tures are achieved with significantly less heat input, and result in 
substantially lower ion temperatures in the topside ionosphere. 
We suggest that this mechanism may be operative under the lim- 
ited low density, refilling conditions in which high ion tempera- 
tures are observed. 

Introduction 


SO 

variation of electron temperature over this range. Predicated on 
the argument that much more thermal energy was transported by 
electrons than by ions, the ion thermal conductivity was not 
treated explicitly. More recently, as numerical simulations have 
become more sophisticated and more numerous, particularly for 
high latitudes, thermal conductivities have again been subjected 
to modification, frequently associated as much with the simula- 
tion conditions as with the plasma physical conditions [e.g. 
Mitchell et al., 1992]. 

The purpose of this report is to propose and test a mecha- 
nism which would decrease ion thermal conduction for a limited 
range of physical conditions. This should increase plasmaspheric 
ion temperatures with less heat input in the low density plasma 
flux tubes where high ion temperatures are most frequently ob- 
served [Comfort et al., 1988). In the next section we describe the 
physical origin of the proposed reduction in thermal conductivity 
and then use this as a basis for formulating a simple mathemati- 
cal model. In the following section, this model is applied in a 
simulation and results with and without the modified thermal 
conductivity are compared. These results and possible refine- 
ments are then discussed and conclusions are summarized. 


In a recent case study examining the nature of heating 
mechanisms required to produce high ion temperatures fre- 
quently observed in low density plasma flux tubes at high alti- 
tudes near the plasmapause. Comfort et al. [1995] found that that 
these temperatures could be readily achieved, provided the heat- 
ing mechanism operated directly on the ions (as described be- 
low), rather than through the thermal electrons. However, the re- 
sulting simulations produced some effects which are not ob- 
served. In particular, ion temperatures in the topside ionosphere 
are up to 5000 K higher than observed, and heavy ion (O'*') con- 
centrations at high altitudes are also much higher than observed. 
While these results are extreme examples, similar problems have 
been found in other case studies, e.g. Horwitz et al. [1990], Cra- 
ven et al. [1995]. Comfort et al. [1995] suggest that these in- 
consistencies between model results and observations could be 
diminished if the thermal conductivity were not so large, so that 
less thermal energy is transported to the topside ionosphere, 
where O"*" is the dominant ion. 

In the early years of in situ observations, modifications to 
the electron thermal conductivity were suggested to explain the 
high electron temperatures observed [Mayr and Volland. 1967; 
Mayr et al., 1973]. These modifications were based on the vari- 
ability of the mean free path over it’s range, and a non-linear 

'Center for Space Plasma and Aeronomic Research 
University of Alabama in Huntsville.Huntsville, AL 35899 
^Space Sciences Laboratory, NASA/Marshall Space Flight Center 
Huntsville, AL 35816 

Copyright 1995 by the American Geophysical Union. 

Paper number 95GL02408 
(X)94-8534/95/95GL-02408$03.00 


Physical Mechanism and Mathematical Model 

The basic idea is most easily understood in the context of a 
mean free path formulation of thermal conductivity. The modifi- 
cation proposed can then be applied to the more complete and 
accurate expression derived from kinetic theory. In the mean 
free path approach, information regarding the thermal properties 
of a fluid for a given location is carried to another location over 
distances of the order of a mean free path (e.g., see derivations in 
Chapman and Cowling [1970. Chapter 6[). These mean free 
paths are limited by particle collisions; in the plasmasphere, these 
are ion-ion (Coulomb) collisions. If densities are sufficiently low 
that collisions are rare, the particle motion will be essentially 
collisionless, and the magnetic moment will be conserved. 

Suppose that the free path is quite long, on the order of the 
length of a flux tube or greater. Those particles which are in the 
atmospheric loss cone will travel to the topside ionosphere and 
lose their energy in collisions with ambient particles. However, 
the many particles which are not in the loss cone will be mag- 
netically mirrored before they undergo a collision. Conse- 
quently, they are returned through the region which they have 
just passed without collision and may well suffer no collision on 
the return trip. As a result, the energy carried away from the lo- 
cation of interest is returned to that location; and, to the extent 
that it was included in the heat to be deposited in other locations, 
the thermal conductivity is in error. The proposal here is 
straightforward: we suggest that only those particles which are 
in the loss cone contribute to the thermal conductivity. This is 
indicated schematically in Figure 1 . While this situation does not 
apply generally, there will be some threshold density, below 
which it will. For the presenf we assume that it holds and inves- 
tigate the consequences. 
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Figure 1. Schematic representation of thermal conduction lim- 
ited to the atmospheric loss cone. Ions which are reflected are 
assumed to transport no heat from the plasmasphere. 

In order to formulate a mathematically simple model, we 
presume that densities are sufficiently low, that ions are essen- 
tially collisionless until they descend to some loss altitude 
helow which they are assumed lost through collisions. The 
magnitude of the magnetic field at that altitude, defines the 
loss cone of the particle distribution at any higher location 
through the conservation of the magnetic moment. That is, par- 
ticles whose mirror field satisfies the relation > fl, will be 
lost. Correspondingly, particles with pitch angles satisfying a < 
ttj, will be lost; and only these particles will contribute to the 
thermal conductivity. We assume an isotropic Maxwellian dis- 
tribution (here the implicit assumption is that on time scales of 
interest there are sufficient collisions to keep the loss cones 
filled, but not enough to invalidate the collisionless assumption 
for particles moving down the field line). Then the fraction of 
particles e, which carries thermal energy out of the plasmasphere 
and into the ionosphere is given by the ratio of the loss cone 
solid angle to the solid angle of the down going flux, 2n. 
From integration of angles over the loss cone, we obtain Clj = 
2n(I-cosai) and e, = I-cosaj From the invariance of the mag- 
netic moment, we can define the loss cone at the position (r.k) 
by 

sin a^^(r,k) I 

B(r.X) 

where r is geocentric distance and X is magnetic latitude. With 
this implicit definition of a^, we can write the expression for e^, 
as 

8^(r,>,)=l- 

V 5/. 

The modified thermal conductivity to be used in the simulations 
is then given by k^= , where is the classical Spitzer- 

Harm thermal conductivity [Spitzer and Harm. 1953]. 

By considering a dipole magnetic field, we can see that 
will vary somewhat with L-shell. But more variation for will 
come from the fact that longer flux tubes with higher equatorial 
crossings will have a wider range of values for Bfr.X). In Figure 
2. we plot values of as a function of geocentric distance for 
different L-shells of a dipole field. From this figure, two aspects 
are readily apparent. First, the L-shell differences are relatively 
minor for a given altitude (providing the L-shell extends to the 
altitude of interest). Second, the modification factor decreases 
rapidly with increasing altitude. It will be noted that this factor 
varies continuously along the field line, being smallest at the 


highest altitudes (equatorial region), with a range of more than 
two orders of magnitude for L > 2. 

Comparisons of Simulations 

In order to determine the realistic consequences for plasma 
thermal structure of using the modified coefficient, we have car- 
ried out simulations with k^., which can he compared with simi- 
lar simulations in which was used. The simulation code is 
the Field Line Interhemispheric Plasma (FLIP) program. This is 
a I-D hydrodynamic code, containing full ionospheric chemistry 
and dynamics; it is described in detail by Richards et at. [1994]. 
To cany out the simulations, we have treated H"*" and O’*" as ma- 
jor ions and He'*' as a minor ion. For the simulations with k^, we 
compute Zj at each step in the calculation and apply it to as 
indicated above. For these simulations we took the loss altitude 
to be 500 km. 

Since the simulations of Comfort et al. [1995] motivated 
this study, we used the same conditions (L = 3, FI 0.7 =231, LT = 
0730, Ap = 3 1 ) so that we could compare results from this study 
with those simulations. In order to produce observed ion tem- 
peratures, using the classical thermal conductivity (as given 
by Banks and Kockarts. Chapter 22), the number (y) of heating 
units ( 1 0’ eV cm'^ s) required to be added directly to the ions was 
8. In those simulations, this heat was distributed evenly over the 
flux tube above 3000 km altitude. (While no specific heating 
mechanism was invoked, it was suggested that Coulomb colli- 
sions with energetic particles or wave-particle interactions might 
cause such heating.) The ion temperature profile corresponding 
to this heat input is shown in Figure 3, labeled as y = 8. The 
other curves are for ic^., with y taking on various values from 8 
down to I. We see that simulated temperatures at the top of the 
flux tube are near those observed when y = 8 using the classical 
conductivity and when y = 1.4 using the modified thermal con- 
ductivity (k^,). Using y = 8 with ic^. results in ion temperatures far 
larger than observed. With k^, considerably less heat is required 
to produce a given temperature, or conversely, a given amount of 
heating produces a much larger temperature increase at high alti- 
tudes. 

To gain a clearer picture of these results, in Figure 4 we ex- 
tract the profiles for y = 1.4 and add the corresponding electron 
temperature profiles. Note that the modified thermal conductiv- 
ity has not been applied to the electrons, since our attention is fo- 
cused on the ions. Comparing the two ion temperature profiles, 
we see that the profile for the modified conductivity case has 
lower temperatures at all altitudes between the top and the bot- 
tom, with smaller gradients near the bottom and larger gradients 



Figure 2. Thermal conduction reduction coefficient as a 
function of geocentric altitude for L-shells 2-6. 






Figure 3. Simulated ion temperature profiles resulting from 
adding heat directly to ions for the case of the unmodified ther- 
mal conductivity and for the modified conductivity with 
Y as a parameter. 

near the top. We can see that the correspondence between ob- 
served and simulated temperatures in the topside ionosphere is 
much closer with the modified conductivity, but the simulated 
temperature at lower altitudes (-800 km) still remains somewhat 
large. The net consequence of all these factors is that less heat is 
conducted out of the plasmasphere to the topside ionosphere. 

The other inconsistency between observations and simula- 
tion results which we were attempting to resolve was the effect 
the ion heating had on the plasma composition. Composition ef- 
fects are seen in Figures 5 a,b. The top panel displays the effects 
of the heating using the classical thermal conductivity, while the 
lower panel shows the effect of the modification. Clearly the 
most dramatic effects are at the top of the flux tube. H'*' densities 
are virtually unaffected, regardless of which k is used, and effects 
on He'*' are also small. However, O'*’ densities have dropped 
significantly. Overall, heavy ion densities are still somewhat 
larger with than observed, but much closer to the observations 
than for the unmodified case. This appears to be a consequence 
of the lower ion temperatures in the topside ionosphere where O'*' 
is the dominant ion species. 

Discussion and Conclusions 

To resolve inconsistencies between simulated ion heating 
results and observations in low density flux tubes, we have pro- 
posed a mechanism which reduces the thermal coupling of the 
plasmasphere to the heat sink of the ionosphere. We suggest that 
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Figure 5. Composition results for the same simulations as in 
Figure 4; observations were made by DEI/RIMS. Panel a shows 
the initial ((7 = 0) H"*", He"*" and 0"*' density profiles and the final 
{7 = 8) density profiles for the unmodified ion thermal conduc- 
tivity (k^//). Panel b shows the heated (7 =1.4) density profiles 
for the modified case (k^). 

magnetic mirroring of quasi-collisionless particles, prevents the 
classical Spitzer-Harm thermal conductivity from providing a 
complete description of the heat flow for sufficiently low densi- 
ties. One question which may be raised is how well the ions in 
our simulation satisfy the quasi-collisionless condition we postu- 
late. To provide insight into this question, we have computed the 
mean time between collisions for protons with energies from 0.5 
eV to 2.0 eV for densities from 50 cm"’ to 500 cm'^, reasonable 
ranges for the outer plasmasphere; these are given in Table 1. 



Temperature (°K) 

Figure 4. Simulated electron and ion temperature profiles for 
heating rates which provide ion temperatures close to those ob- 
served at high altitudes for cases with (7 =1.4) and without (7 = 
8) modified ion thermal conductivity. High altitude observation 
is from DEI /RIMS and low altitude observation is from 
DE2/RPA [from Horwitz et ai, 1990]. 


We have also estimated transport times for protons with these 
same energies by considering the limiting case of particles mov- 
ing along the field line for L-shells from 2 to 5, from initial posi- 
tions at the equator and at latitude 45° down to the collision alti- 
tude (500 km); no forces are included in these simple estimates. 
Table 2 shows the resulting transport times. For the conditions 
of the simulation above (n = 100 cm'\ L = 3), the ratio of the 
collision time to the transport time is about 2 for protons initially 
at the equator, which is in the appropriate range for the suggested 
model. 

Table 1. Proton collision time (sec) as a 

function of energy and density 


Energy (eV) 

n (cm'^) 

0.5 

1.0 

2.0 

50 

2178 

6161 

17430 

100 

1089 

3081 

8713 

200 

545 

1540 

4357 

500 

218 

616 

1743 
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Table 2. Proton transit time (sec) from initial point (L,A.) 
to 500 km altitude as a function of proton energy 


Energy (eV) 


L 

0.5 

1.0 

2.0 

(>^) 

(0°) 

(45°) 

(0°) 

(45°) 

(0°) (45°) 

2 

1117 

* 

790 

* 

558 * 

3 

2029 

354 

1435 

250 

1015 177 

4 

2934 

700 

2074 

495 

1467 350 

5 

3836 

1044 

2712 

738 

1918 522 


This has basically been a ‘proof of concept’ investigation to 
determine whether or not the proposed mechanism would pro- 
duce the type of effects the observations require. As such, it 
applies to physically limited conditions in which the densities are 
such that the plasma is in a quasi-collisionless to collisional re- 
gime. The simulations suggest that if the mechanism were op- 
erating, it would have the qualitative characteristics necessary to 
bring simulation results closer to observations. However, a more 
refined quantitative study (which we plan on carrying out) 
should consider additional aspects which are quantitatively im- 
portant. For example, as seen in Figure 5 above, the number 
density varies along the field line so that the mean free path also 
varies with position. Similarly, the longer path lengths for the 
higher L-shells, make satisfaction of the collisionless criterion 
more restrictive for a given particle population. In addition, the 
particle path length is actually longer than the field line distance 
(guiding center path length), depending on the particle’s equato- 
rial pitch angle; the field line length is a lower limit. These ef- 
fects should be folded into a density dependent modification to 
the thermal conductivity. A density dependence was developed 
by Mayr and Volland [1967] for electrons, but it had no relation 
to the effect of mirroring described in this report. 

Since the Coulomb cross section is energy dependent, par- 
ticles with different energies will have different mean free paths, 
as indicated in Table 1. For the quasi-Maxwellian velocity dis- 
tributions expected for these conditions, it is possible that the 
high energy part of the distribution could satisfy the collisionless 
criterion, while the low energy part would not. A kinetic treat- 
ment will be required to examine this aspect in detail. 

It appears likely that the mirroring mechanism proposed here 
can play a significant role in reducing the thermal coupling of the 
plasmasphere to the ionospheric heat sink, for the restricted 
conditions indicated. While the simple model used to test the 
idea is applicable under only these conditions, the physical de- 
coupling of the ionosphere from the plasmasphere for particles 
outside the loss cone, due to magnetic mirroring, will occur for a 
broader range of lower densities. This mechanism is basically a 
kinetic phenomenon, and is beyond the means of even higher or- 
der hydrodynamic models to include. In this respect, the simple 
model suggested can be considered an attempt to extend the 
range of the hydrodynamic models into the quasi-collisionless 
regime with a physically-based, ad hoc modification. More de- 
tailed examination, with the aid of case-study observations and a 


kinetic simulation code, are needed to adequately determine the 
appropriate manner in which to employ this mechanism in a hy- 
drodynamic code. 
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